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Abstract 

A generalization of the procedure to study shapes and stability of algebraic nuclear models 
introduced by Gilmore is presented. One calculates the expectation value of the hamiltonian 
with respect to the coherent states of the algebraic structure of the system. Then equilibrium 
configurations of the resulting energy surface, which depends in general on state variables 
and a set of parameters, are classified through the Catastrophe theory. For one and two- 
body interactions in the hamiltonian of the interacting boson model-1, the critical points 
are organized through the Cusp catastrophe. As an example, we apply this Scparatrix to 
describe the energy surfaces associated to the Rutcnium and Samarium isotopes. 


1 Introduction 

The geometry of algebraic nuclear models can be studied by means of the time-dependent varia- 
tional principle [1], [2]. This formalism provides us with a classical limit of the nuclear model, in 
particular we are mainly concerned with the static properties of the hamiltonian function (energy 
surface) associated to the considered algebraic nuclear model. In general these hamiltonian func- 
tions depend on state variables and a set of parameters, then the appropriate mathematical tool 
to determine the most general behaviour of their equilibrium configurations is the Catastrophe 
formalism [3]. 

A connection between the interacting boson modcl-l(IBM-l) [4] and the geometrical approach 
of Bohr-Mottelson [5] was done by expressing the IBM-1 hamiltonian in terms of shape variables. 
This can be achieved by means of the intrinsic boson states defined by [6] or by the corresponding 
coherent states [2]. Analysis of shape and phase transitions in this model have been done by [7], 
[8]. In this work we apply the procedure introduced in Ref. [2] to the interacting boson model, 
but for the general hamiltonian of one and two-body central interactions involving s and d bosons 
[4] and determining its associated Scparatrix. We show that the equilibrium configurations can be 
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classified through the Cusp catastrophe, this means that two parameters are enough to describe 
the most general energy surface. Therefore this analysis generalize thpse presented previously 
[4], in which only transitions between pairs of exact SU( 5), 0(6) and SU (3) symmetries are 
considered. In the last decade, effective hamiltonians of the IBM-1 have been used to describe 
energy spectra and transition probabilities of chains of isotopes and isotones [9], [10]. In particular, 
the effective hamiltonians for Ru [11] and Sm [12] isotopes were determined, i.e., the best choice of 
the parameters of the general IBM-1 hamiltonian that reproduced the corresponding experimental 
data. Using these effective hamiltonians we construct their energy surfaces and show that their 
critical points follow a curve in the parameters space organized by the Cusp Separatrix. This let 
us to know: i) How many equilibrium configurations yield the system and ii) If the behaviour of 
the model around the critical points may or may not be approximated by an harmonic oscillator. 

In the Second Section we review how the energy surfaces can be determined considering a 
hamiltonian constructed in terms of the generators of a Heisenberg- Weyl algebra. In the Third 
Section a brief summary of the IBM-1 is presented. In the Fourth Section, an analysis of the 
shape and stability of the most general energy surfaces of the IBM-1 is made, also the curves 
associated to the Ru and Sin isotopes are plotted in the parameters space. Finally some remarks 
are indicated in the last section. 


2 Energy Surfaces of Algebraic Models 

The energy surfaces (ES) of algebraic models can be determined by means of the coherent states 
of the associated algebraic structure of the hamiltonian. As an example, a hamiltonian written in 
terms of the generators of a Heisenberg- Weyl algebra is considered, i.e., 

H = cq b + C 2 ^ b^) + Cj (6^ + b) , (1) 

where the operators 6* and b satisfy standard creation and annihilation commutation relations. 
Although this hamiltonian can be solved analytically by means of a Bogoliubov transformation, 
we use it to illustrate the procedure to construct the coherent states and the ES of an algebraic 
model. 

The coherent state is defined by the action of the raising generator on the vacuum state [1] 

|q) = exp(a* 6*)|0) - (2) 


The Baker-Campbell-Hausdorff formulas can be used to calculate the overlap of two coherent 
states and the coherent state representation of the creation and annihilation boson operators 


(a|a) = exp (a' ot*) , 
(a\b=-^-(a\, (a|6 t = a(a| 


Then the energy surface is given by 


(aM/fda) 

E(a) = Inn Vt\ • 

a —¥ ot (ft | fV/ 


(3) 

(4) 

(5) 
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Substituting the expressions (1), (3) and (4) into the previous one, one gets the result 

E{p) = (c 0 + 2 c 2 )p 2 + 2ci/> , (6) 

where the reality condition a = a* = p was also used. 


3 The Interacting Boson Model 


In 1975 this model [4] was introduced to describe collective properties of even-even nuclei far from 
closed shells through the interactions between two kinds of bosons, one with angular momentum 
L = 0 (the s-boson) and another with angular momentum L = 2 (the d-boson). The six possible 
boson states give rise to a U( 6) group structure. The bosons represent pairs of fermions, the 
s-boson reflects the strong pairing attraction of identical nucleons whereas the d-boson is a result 
of the weaker J — 2 + attraction [4]. Therefore nuclei are pictured as systems of s and d bosons, 
whose number is equal to half the number of the valence nucleons, the core being considered inert. 
When a shell is more than half full, hole-pairs are counted instead of particle-pairs. 

The most general one and two body hamiltonian that conserves the total number of bosons is 


Hibm = £3 N t + £4 N,i 4 - £ ^v'2ZTl[[d t xdtjW x[dxd]M]M 

0,2,4 Z 

+ 4= ( [(d f X d*] l 2 ) X d] t°l S + S* [d f X [d X d] M] t°l ) 

v2 

+ | (((d l xd']l 01 » 2 + » ,2 [dxd|i°l) 

+ v^5u 0 s’ s [d* x d] + «2 s t2 s 2 , (7} 


where the sets of boson operators s, s* and dj t , d fl satisfy the following, different from zero, com- 
mutators 



( 8 ) 


Now, we construct the coherent state of a six dimensional harmonic oscillator, following the 
procedure indicated in the previous section. However in this case the associated group is compact 
and then we restrict the exponential to only one term of the Taylor series expansion, 


|W, o) = • 

A* 


( 9 ) 


where is the normalization constant. Evaluating the corresponding Eq.(5) one arrives to the 
formulae for the energy surface of the model [8], t.e., 

E{(3 , 7) = N € ^ ff 2 ) + + ( a i/? 4 + «2/? 3 cos37 + a 3 /? 2 + u 2 ) , (10) 

where it was used that the laboratory variables can be expressed in terms of two intrinsic 
parameters (3 and 7 plus three Euler angles. Besides as the energy surface is a rotational invariant 
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all the dependence in the Euler angles disappear. The parameters aj, 02, a 3 , and e of the Eq.(10) 
are combinations of those that appear in the IBM-1 hamiltonian (7) 


ai = 

c 0 c 2 9c 4 
10 + 7 + 35 ’ 

(11) 

a 2 = 

~{ts V o ' 

(12) 

a 3 = 


(13) 

€ = 

U ~ G • 

(14) 


4 Shapes and Stability of Energy Surfaces 

The energy surfaces define functions of state variables and a set of parameters, and the Catastrophe 
theory is used to analize their equilibrium configurations. This formalism let us organize all the 
possible shapes of the ES into well defined separated regions of the parameters space. 

To illustrate how this is done, we consider the potential energy surfaces (PES) of the simplest 
version of the Generalized Collective Model [5], i.e., 

V(P, 7; C 2 , C 3 ) = P 4 - C 3 P 3 cos3j + C 2 p 2 . (15) 


The equilibrium or critical points are determined by solving the equation y V" (/?, 7) = 0. The 
results arc given by (0,0) and (/?oi7o)> with 


A> = 


3C 3 ± yj$Cl - 32C 2 
8 


7o = 0, ?r/3 . 


(16) 


The set of degenerated critical points defines a locus in the parameters space which is called the 
Scparatrix. This can be obtained through the determinant of the Hessian matrix or by other 
procedures, in this case it is immediate that the critical points are double degenerated if the 
parameters satisfy the expressions 


9 C 3 — 32 C 2 = 0, (17) 

C 2 = 0 , C 3 7^ 0 . (18) 

For the expression (17) the critical points arc localized in Pq = 3C 3 /8 while for (18) in /? 0 = 0. 
Besides, it is straight to prove that if C 2 = 0 and C 3 = 0 the critical point is triple degenerated 
and localized at p 0 = 0. The Scparatrix of the system is shown in Fig.l. It divides the space in 
regions each characterizing a typical shape yield by the model. By means of the transformation 
P = y + C 3 /4 the Scparatrix is taken to the canonical form of the Cusp catastrophe. 

Now we study the equilibrium configurations of the energy surface associated to the IBM-1, 
which is given in Eq.(10). Then we calculate the critical points by taking the derivatives with 
respect to P and 7 variables. It is straightforward to sec that the critical points correspond only 
to 7 = 0 (prolate case) or 7 = tt/ 3 (oblate case). Therefore we can restrict to the prolate case, 
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without losing generality. The energy surface (10), with 7 — 0, can be rc-wiitten in terms of the 
following parameters 


ri 


r 2 


2 u 2 (N - 1 )-(N- l)fl 3 - 6 
2 a\(N - 1 ) — (N - 1)«3 + « 
2 o 2 (iV - 1 ) 

2ai(7V — 1 ) — (N — 1)03 + c 


( 19 ) 

(20) 


and it takes the form 

m = (i + pjt ^ + r,02{l>2 +2) ~ r2/j3 l • (21) 

where E(P) = E{P , 7 = 0)/N - e 3 - u 2 (N - 1). One has to notice that the oblate ease can be 
regained by interchanging r 2 by —r 2 or equivalently a 2 by —a 2 . 



FIG. 1. Separatrix for the Bohr Mottelson Hamiltonian 

To find the extrema in /? of the Eq.(21) one needs to solve the equation 

P ( r 2 /? 3 + 4 /? 2 - 3j - 2 ^ + 4r 1 ) = 0 . ( 22 ) 

From this expression, one determines the locus in the parameters space ( v 2 , 77 ) of degenerated 
critical points. Then the Separatrix of the model is defined by the curves 


77 


»'! 


(16 4 - 9» 2 ) 3 ^ 2 32 

± 54 ?-| 27 j -| 

0 . 


( 23 ) 

( 24 ) 


This Separatrix is sliovVn in Fig.2 and it corresponds to the Cusp catastrophe^ although it does not 
has the canonical form. 
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FIG. 2. Scparatrix for the IBM-1 Hamiltonian 

Now we applied the results to find the shapes and stability of the Ru and Sm isotopes. For 
the Ru case, one has the number of valence protons pairs, N , = 3. As we consider isotopes with a 
mass number varying from A = 98 — 110, the corresponding number of neutron bosons runs from 
N„ = 2 to N v — 8. Thus the size of the space in the IBM-1 is determined by the total number 
of bosons, which is the sum of the numbers of proton and neutron bosons. In the Sm case, one 
has the number of proton bosons, N r = 6, and as we take into account A = 148 — 160, the 
number of neutron bosons runs from N„ = 2 to N u = 8. An analysis of their energy spectra and 
electromagnetic transitions using the hamiltonian (7) is made in [11] and [12]. For both isotope 
chains, the parameters used are presented in the first eight columns of Table I. Substituting these 
parameters into the equations (11) to (13) we get the values of ai ,a 2 and a 3 . These are given in 
the last three columns of Table I, by means of which one can easily construct the corresponding 
energy surface of each nucleus. 

To find the region of the Separatrix , Fig. 2, where the different isotopes are localized one 
calculates the parameters and r 2 through the equations (19) and (20), as functions of the total 
number of bosons. For the Ru , one gets the expressions 


_ 990.2+ 146.2(JV- 1) 
n ” 990.2 + 40.2(N — 1) ’ 

while for the Sm isotopes the parameters are 

_ 2171.2 + 258.3 (N — l) 
n “ 2171.2 + 151.2 (N — 1) 


>2 = 0 , 


86 (N — 1) 
2175— 151 (N - 1) ' 


The localization of the joints (25) and (26) are shown in Fig. 3 and Fig 4, respectively. 


TABLE Parameters, in KeV_ units, used to describe the Ru and Sm isotope chains 
t Co c 2 C 4 u 0 «2 »0 V2 Ol | 02 I oT" 

Ru 990.2 -185.4 -77.4 -.0.4 -53.0 23.3 0 -52.1 -29.7 ~0 ^52 1 


Co 

C 2 

C 4 

Wo 

«2 

vo 

t>2 

a i 

a 2 

fl 3 

-185.4 

-77.4 

-.0.4 

-53.0 

23.3 

0 

-52.1 

-29.7 

0 

-52.1 

-613.7 

-318.8 

-377.6 

227.4 

0.4 

89.9 

-33.0 

-204.1 

-43.0 

-256.9 
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FIG. 4. Plot of the points (26) associated to the Sm isotopes chain. 

5 Conclusions 

For the Bohr-Mottclson model the yielded shapes and equilibrium configurations arc classified by 
the Scjmratrix of Fig.l. One can identify four regions*, (i) Above the parabola, the PES have one 
minimum at P = 0. ii) Between the parabola and the C 2 axis, in the PES appeal additionally 
a second minumum at P ^ 0. iii) Below the C 3 axis the PES have two mimima at p / 0 and 
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the maximum occurs at (3 = 0. The PES around these critical points can be approximated by 
an harmonic oscillator, iv) On the locus of points that define the Scparatrix , we have shape 
transitions and in the vecinity of the critical points (16) the PES cannot be approximated by 
quadratic functions. It is important to remark that the PES have a mirror symmetry along the 
C3 axis, which physically represents transitions from prolate to oblate shapes. 

For the IBM-1 hamiltonian one gets (see Fig. 2): (i) For positive values of r { one has two 
regions, above the curve the ES present one minimum at /? = 0 and below the curve they have 
a second minumum at (5 ^ 0. ii) For negative values of r x one has again two regions, above the 
curve the ES are built with two minima at /? / 0 and a maximum at (3 = 0 and below of it, 
the ES have one mimimum point at (3 ^ 0 and a maximum at /? = 0. For the critical points 
mentioned above the energy surfaces can be approximated by an harmonic oscillator, iii) Finally 
on the Scparatrix, there are shape transitions, and in the vecinity of their critical points the ES 
cannot be approximated by quadratic functions. It is important to remark that these ES have 
also a mirror symmetry along the r 2 axis, representing transitions from prolate to oblate shapes. 
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